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P-representability is a necessary and sufficient condition for separability of bipartite Gaussian 
states only for the special subset of states whose covariance matrix are Sp(2, R) ® Sp(2, R) locally 
invariant. Although this special class of states can be reached by a convenient Sp(2, R) ® Sp{2, R) 
transformation over an arbitrary covariance matrix, it represents a loss of generality, avoiding infer- 
ence of many general aspects of separability of bipartite Gaussian states. 
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I. INTRODUCTION 

In the recent years the question to whether a given 
quantum state is separable or entangled has become cen- 
tral to the quantum information and to the quantum op- 
tics communities. Mostly because fault-tolerant quantum 
information protocols, such as quantum computation and 
quantum teleportation are completely dependent on the 
ability to prepare pure (or close to pure) entangled states 
[1]. Recent attention however has been centered on 
continuous variable versions of quantum communication 
protocols, such as the unconditional quantum teleporta- 
tion [2,3], whose efficiency rests on the ability to gener- 
ate entangled states of systems with infinite dimensional 
Hilbert spaces. Bipartite systems with finite Hilbcrt 
space have been exhaustively investigated in order to 
achieve a precise quantification of entanglement. Peres 
[4] and Horodecki [5] demonstrated that a necessary and 
sufficient condition for separability of bipartite systems 
with Hilbert spaces of dimension < 2 ® 3 is the positivity 
of the partial transpose of the system density matrix. On 
the other hand, algebraic similarities between bipartite 
states with Hilbert space of dimension 2 (g> 2 and bipartite 
Gaussian states (described by 4 x 4 covariance matrices) 
allows the extension of the positivity criterion to those 
special continuous variable states as firstly developed in 
Refs. [6,7], and considered afterwards in discussions on 
entanglement of Gaussian bipartite states (e.g. [8-18]). 

Of particular importance is the connection between 
Glauber P-representability of a bipartite quantum state 
and separability [10]. A P-representable state is the one 
that is represented by a positive Glauber P distribution 
function P(a,P), which is less(or equally)-singular than 
the delta distribution, such as 

p = J da 2 df3 2 P{a,(3)\a,(3)(a,f}\. (1) 

Under this condition P(a, (3) assumes the structure of a 
legitimate probability distribution function over an en- 
semble of states, allowing the connection between sep- 



arability and classicality. However, although any P- 
representable bipartite state is separable, as can be im- 
mediately seen by the P representation definition, the 
inverse is not necessarily true. P-representability and 
separability are completely equivalent only for Gaus- 
sian states with locally Sp(2, R) ® Sp(2, R) invariant 
covariance matrices. Since any covariance matrix can 
be brought to this invariant form under appropriate 
Sp(2,R) (g) Sp(2,R) transform, P-representability and 
separability have been misleadingly accepted as one-to- 
one equivalent properties of bipartite Gaussian states. 
The purpose of the present paper is to give a complete 
classification of the set of all bipartite Gaussian separable 
states (BGSS). Particularly we show that P-representable 
Gaussian bipartite states form a subset of BGSS with lo- 
cally Sp(2, R) ® Sp(2, R) invariant form. In Section II we 
begin by a revision of some necessary properties of bipar- 
tite Gaussian states and in Sec. Ill we give the necessary 
and sufficient conditions for the state to be separable. In 
Sec. IV we discuss the P-representability of those states 
and show that they actually form a subset of the sep- 
arable states. In Sec. V we describe the the unitary 
Sp(2, R) (g> Sp(2, R) map that connects the two sets and 
finally in Sec. VI a conclusion encloses the paper. 

II. BIPARTITE GAUSSIAN STATES 

Any bipartite quantum state p is Gaussian (see e.g. 
[17,19]) if its symmetric characteristic function is given 
by 

C(r,) = Tr[D(r,)p] = e -^ +v ", (2) 

where D{rf) = e _T,tEv is a displacement operator in the 
parameter four-vector 77-space, with 

»7 T = (ViiViiVhrn) , v + = (a{, 01,4,02) , (3) 

and 
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where a\ (a\) and a 2 (<4) are annihilation (creation) op- 
erators for party 1 and 2, respectively. V is the Her- 
mitian 4x4 covariance matrix with elements Vij — 
(— iy +: >({vi, Vj})/2, which can be decomposed in four 
block 2x2 matrices, 



Vi C 

ct v 2 



(5) 



where Vi and V 2 are Hermitian matrices containing only 
local elements while C is the correlation between the two 
parties. Any covariance matrix must be positive semidef- 
inite (V > 0), furthermore the generalized uncertainty 
principle 



we can further simplify the generalized uncertainty (9) 
and (10) to 



ni > \ \mi\ 2 + -, 



(12) 



and 




to. 



+ |to 2 -c| 2 , (13) 



respectively, with 
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nt (|to c | 2 + |to s | 2 ) — m c m s m\ — to*to*toi, (14) 
2nim*m c — to 2 to* — (to*) mi, (15) 



h ~ 4 - \ m i\ 



(16) 



V + -E>0, 



(6) 



must also be applied. Those general positivity criteria 
can be decomposed into block using matrix positivity 
properties. There are many ways to check a matrix pos- 
itivity, such as the positivity of the matrix determinant. 
However this is not a necessary condition. A reliable and 
convenient way to check the covariance matrix positiv- 
ity is through the following block Schur decomposition 
[20]: Any Hermitian matrix is positive if and only if any 
principal block matrix is also positive, or, if its upper left 
block and the block's Schur complement are also positive. 
So that for the covariance matrix (5), V > only if 



Vi >0, 

and the Schur complement of Vi 

S(Vi) = V a - C^Vi^C > 0. 



(7) 



(8) 



It is interesting to observe that Schur complements of 
block matrices representing Gaussian states covariances, 
such as above, embodies a manifestation of a physical op- 
eration when considering partial projections onto Gaus- 
sian states [21,22]. 

Through the Schur decomposition the physical positiv- 
ity criterion (6) applies only if 



Vi + -Z > 0, 



(9) 



and 



(V 2 + ±Z) - Ct (Vx + iz) C > 0. (10) 



By explicitly writing 
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(11) 



III. SEPARABILITY BOUNDS 



By mapping the positivity necessary and sufficient con- 
dition for dimension 2 ® 2 to bipartite systems of infi- 
nite dimension, Simon [6] has discovered an elegant ge- 
ometrical interpretation of separability in terms of the 
Wigner distribution function for the density operator. 
The Peres-Horodecki separability criterion in the Simon 
framework reads: if a bipartite density operator is sepa- 
rable, then its Wigner distribution necessarily goes over 
into a Wigner distribution under a phase space mirror 
reflection. The separability criterion can be understood 
as a valid Wigner-class-conservative quantum map un- 
der local time reversal. Following [6] a necessary and 
sufficient condition for a Gaussian quantum state to be 
separable, 



(17) 



is that its covariance matrix must satisfy 



V+iE>0, 



(18) 



under a partial phase space mirror reflection (partial Her- 
mitian conjugation) V = TVT : 



Tv = v T = 



/ai\ 



(19) 



with 



T = 



I 

X 



X = 
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(20) 
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otherwise the state is entangled. Similarly to the gen- 
eralized uncertainty (6) decomposition, the separability 
condition (18) is satisfied if and only if (9) and 

^XV 2 X+^Z^ -XC 1 (Vi + iz^j CX>0, (21) 
are both satisfied, which explicitly implies in (12) and 




d 



+ 1 



m 2 



(22) 



respectively. We call the set of states p that fall inside 
(12) and (22) the set BGSS of all bipartite Gaussian sep- 
arable states S. Any state that does not follow those 
inequalities is entangled, being it pure or not. Remark 
also that purity is only reached when the equalities in 
(12) and (13) hold. 



IV. P-REPRESENTABILITY OF GAUSSIAN 
STATES 

The very definition of a separable state (17) can be 
written in a coherent state representation through the 
Glauber P-function (1), but it is not obvious that P(a, (3) 
is a legitimate probability distribution function. That is 
only reached if the state is P-representable, i.e., if the P- 
function is non-negative and less (or equal) singular than 
the delta distribution. In terms of the covariance matrix, 
a quantum state is P-representable [19] if 



v-ii>o, 



(23) 



which in terms of the upper left block matrix and its 
Schur complement writes as 



Vl -il>0, 



and 



Thus 



V 2 -I 



Vi 



C> 0. 



(24) 



(25) 



States that follow (26) and (27) form the set of all bipar- 
tite P-representable Gaussian states P. 

From an operator formalism for the density matrix, 
Englcrt and Wodkiewicz [10] have recently stated that 
P-representability is equivalent to the separability con- 
dition, for the specific symmetric situation where mi = 
to 2 = m s = 0, rii = n 2 = n, and m c = to, which indeed 
set P ^ S as we see bellow. The generality of their state- 
ment is justified only if Sp(2, R) ® Sp(2, R) local opera- 
tions are used to bring those parameters to the special 
symmetric class described above (see also [19]). How- 
ever, this particular situation does not represents a total 
equivalence between § and the set of all P-representable 
states. In general the P-representability conditions, (26) 
and (27), are more restrictive than the separability ones, 

(12) and (22), respectively, as we now investigate. 

Firstly observe that (12) is less restrictive than (26), 
equaling only for \m\\ — or \m\\ — ► oo, being enough 
to check if (27) dominates over (22) for the simplest 
| mi | = situation. For that we make use of the knowl- 
edge that (13) is always stronger than (27), including 
the situation where d = 0, i.e., m = 1/2. In such a case, 
the comparisons of the (27) lower bound to (22) and to 

(13) are equivalent and thus if (22) is violated so is (13). 
These inequalities must satisfy 



(|m 2 | + |m c | 2 )(|m 2 | + |m s | 2 ) >0, 



(31) 



and since the quantities involved are always strictly pos- 
itive the criterion (31) is always satisfied. The equal- 
ity however occurs only if |m 2 | = |m c | 2 = or |m 2 | = 
|m s | 2 = 0, which then set the equivalence P ^ § for the 
two following special Sp(2, R) <g> Sp(2, R) invariant forms 
forV: 

Invariant form 1: 



(32) 
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(26) 



Invariant form 2: 
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(33) 



Special forms 1 and 2 arc locally Sp(2, R) Sp(2, R) in- 
variant covariance matrices that form the P ^ S subset. 
The separability and thus P-representability criterion is 
then reduced to 
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U1 ~ \) V 2 ~ \ } - ' TOC ' 2 



for the special form 1, and to 



1\ ( 1 



ni ~ 2 ) \ 2 ~ 2 ' - ' ms ' ' 



(34) 



(35) 



for the special form 2. The physical condition of exis- 
tence of a general bipartite Gaussian state of the form 1 
or 2 writes 



- 2) \ n2 + \ ) - ' mc ' 2 



or 











f n 2 H 









(36) 



(37) 



respectively. 




Fig.l. Typical Separability (S) and P-representability 
(P) boundaries for mi = 0.5 and mj = 1. The shaded 
area where the P-fold is lower than the S-fold does not 
represent physical quantum states. 

Remark 1: There are P-representable Gaussian oper- 
ators that violate (13), which however do not represent 
any valid positive definite quantum state. As an exam- 
ple in Fig. 1 through comparison of the limiting bounds 
(22) and (27), assuming all real coefficients and setting 
nil = 0.5, and m 2 = 1. Only those states that lay in or 
above the separable class are valid P-representable sepa- 
rable states. 

Remark 2: The special symmetric situation depicted 
in Ref. [10,13] for the two-mode thermal squeezed state, 
where mi = rri2 = m s = 0, n\ = n 2 = n, and m c = m, 
is a particular example of the specific form 1, and thus a 
separable state in this case is always P-representable. 



V. § <=> P MAPPING 

Any general covariance matrix can be mapped into one 
of those invariant forms under appropriate Sp(2, R) £§> 
Sp(2, R) transform. In other words, it is possible to map 
S into P such that 



Ps P = UlpqUl 1 , 



(38) 



be the state obtained by the local unitary transform 
Ul = U\ ® U2 over a general bipartite Gaussian density 
operator p<j assuming 



U L vUl l = S L v, S L = ( q 1 g ° 2 



(39) 



with the condition Sl 1 = ESl t E. The new symmetric 
characteristic function writes 



C Sp (T7) = Tr[D( V )p Sp ] = Tr[U^D( V )U L p c 



with 



V 5p = Sl+VSl. 



(40) 



(41) 



The transformed covariance matrix writes as (5), but 
with new block elements 

V; = S i t V i S i , C' = Si t CS 2 . (42) 

Assuming a local Sp(2, R) transform as 



e^coshfli e i(p > sinhfl, 
e"^ 1 sinh 9; e"^ 1 cosh ft 



(43) 



the condition to bring V to the invariant form 1: 

v\ p c 

vi p* c 

p c V2 

Pc v 2 



(44) 



is obtained by setting 4>i + tfi = —p% + 7r and tanh 2ft = 
I m j I /rii = |m c |/|m s |, for i — 1,2, respectively, where 



-ifii 



TO,/|mj|, and assuming (for \m c \ > |m s |). 



Now the condition to bring V to the invariant form 2: 

(45) 



v\ p s 

vi p* 

p* s u 2 

p s v 2 



is immediately attained if <pi + tp, = —pi + tt also, but 
now with tanh 2ft = \mi\/rii = \m s \/\m c \ (assuming 
\m c \ < |m s |). Since both Vi and Vi are proportional 
to the identity, they do not change under unitary local 
rotations and the two invariant forms are then connected 
through those operations. As such, the last two con- 
ditions on \m c \ and \m s \ can be waved by appropriate 
rotations. 

The new transformed elements are 



p s = e -<(*i-te). 



■\/\m s \ 2 ~ |m c | 2 , 



(46) 
(47) 
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VI. CONCLUSION 



^ = e-^+«f^WK| 2 -KI 2 , (48) 
\m c \ 

(for \m c \ > |m s |), which then turn explicit the four 
invariants of the Sp(2, R) (g) Sp(2, R) group: I\ = 
detVi, I 2 = detV^, h = det C, and I 4 = 
Tr [ViZC'ZV^Z(C') t Z]. 

The general Sp(2, R)®Sp(2, R) transformation (41) of 
the (43) form is reached through the squeezing operation 
U L = U x ® U 2 : 

over the bipartite Gaussian state pc, with \ni\t = 9i = 
2ri, the squeezing parameter associated with the trans- 
formation on the mode i and e l ^ i = m/\m\. An 
important result is that while all the BGSS set can 
be mapped into the P-representable set by suitable 
Sp(2, R)®Sp(2, R) transforms, it is not possible to resti- 
tute the original matrices Vi and V 2 with unitary ro- 
tations. That is only reached applying over the squeez- 
ing operation. This is immediate from the two invariant 
forms. Since both covariances reduced matrices V^ and 
V 2 are proportional to the identity, unitary rotations 
transform the invariant forms among themselves. 

Remark 3: Through the § ^ P mapping, we have 
reached the special subset of locally Sp(2, R) ® Sp(2, R) 
invariant forms. However the general separability con- 
dition can possibly be set equivalent to the special P- 
representable subset under an appropriate nonlocal oper- 
ation forming the mapping 8 — > P. Let us consider again 
the condition for separability (18). It can be equivalently 
written as 

V + ^TET > 0. (50) 

Now let Uml be a nonlocal operation: 

U N Lvtf NL = Mv, (51) 

where M is a general transformation matrix: M G 
Sp(4,R). Such a general M, when acting on (50) must 
leave V invariant in form (V), while M^TETM must 
go necessarily to —I, such that (50) writes as 

V - h > 0, (52) 

which is the transformed P-representable subset condi- 
tion. The Stone - von Neumann theorem provides that 
if M exists it must be unitarily implementable [23]. 
Nonetheless, finding the corresponding Unl operator 
may not be a simple exercise and we leave this point 
for future research. 



In conclusion, we have derived a complete descrip- 
tion of bipartite Gaussian separable states, and have 
proved that P-representable states form a subset of 
the set of all bipartite Gaussian separable states, ex- 
istent only under special symmetry of the covariancc 
matrix. We can state that for positive definite bipar- 
tite Gaussian operators, which describe physical quan- 
tum states, P-representability is a necessary and suffi- 
cient condition for separability only for the subset of lo- 
cally Sp{2 1 R) <g) Sp(2,R) invariant Gaussian states [24]. 
In General PcS. 
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